Transmission and reflection of strongly nonlinear solitary waves at granular interfaces 
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The interaction of a solitary wave front with an interface formed by two strongly-nonlinear 
non-cohesive granular lattices displays rich behaviour, characterized by the breakdown of continuum 
equations of motion in the vicinity of the interface. By treating the solitary wave as a quasiparticle 
with an effective mass, we construct an intuitive (energy and linear momentum conserving) discrete 
model to predict the amplitudes of the transmitted solitary waves generated when an incident 
solitary wave front, parallel to the interface, moves from a denser to a lighter granular hexagonal 
lattice. Our findings are corroborated with simulations. We then successfully extend this model 
to oblique interfaces, where we find that the angle of refraction and reflection of a solitary wave 
follows, below a critical value, an analogue of Snell's law in which the solitary wave speed replaces 
the speed of sound, which is zero in the sonic vacuum. 
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PACS numbers: 45.70.-n, 61.43.Fs, 65.60.+a, 83.80.Fg 



The study of solitary waves has over the years led 
to a paradigmatic shift in our understanding of many 
body phenomena characterized by anharmonic effects 
that manifest themselves in exotic electronic [I[ and me- 
chanical states 0, H| . A concrete and technologically rel- 
evant 0, [B[ arena to study strongly non-linear mechan- 
ical waves is the sonic vacuum 0, 0] ~ a paradigmatic 
example of which is an aggregrate of grains just in con- 
tact. Owing to the vanishing speed of linear sound, even 
the tiniest strains propagate as supersonic solitary waves, 
non-linear periodic waves and shock-like waves depend- 
ing on conditions of loading. However, the differential 
equations describing the propagation of mechanical dis- 
turbances around the state of sonic vacuum are both non- 
linear and generally not integrable, making it difficult to 
model solve them analytically. Moreover, the continuum 
approximation itself may fail in the vicinity of a sharp 
granular interface, where the discrete nature of the gran- 
ular medium dominates. 

So far most studies of strongly non-linear granular in- 
terfaces have concentrated on what happens when a soli- 
tary wave initially propagating in a chain with mass m\ 
reaches an interface where the particles mass suddenly 
changes to mi- Depending on the ratio A = qual- 
itatively different behaviours are observed. When the 
solitary wave moves from a lighter to a denser medium 
(A > 1), most of its energy gets divided into a reflected 
and a transmitted pulse, whose respective amplitudes 
can be estimated using the conservation of linear mo- 
mentum and energy 0j H 0- By contrast, when the 
incident solitary wave moves from a denser to a lighter 
medium (A <C 1), a train of (multiple) solitary waves is 
generated in the lighter medium |2j, [6|. In this case, a 
direct application of the two conservation laws (momen- 
tum and energy) is not sufficient to predict the ampli- 



tude ratios of the multiple solitary waves. Experimental 
studies of granular chains comprised of steel and polytc- 
trafluoroethylene (PTFE) particles (A = 0.27) have shed 
light on the discrete mechanism responsible for the gen- 
eration of the train of solitary waves in the PTFE light 
chain. Most of the collective motion carried by the inci- 
dent solitary wave (propagating in the chain of stainless 
steel particles) is converted into the motion of a single 
interfacial steel particle [1, ||| [ll[ . 

In this Letter, we turn to the hitherto unexplored two 
dimensional problem of determining the reflection and 
transmission of a strongly non-linear solitary wave-front 
incident upon an interface between two hexagonal lat- 
tices both in a sonic vacuum, but with different par- 
ticle masses. We treat the solitary waves as quasipar- 
ticles with an effective mass and model the interaction 
with a two-dimensional granular interface, by assuming 
an energy and linear momentum conserving scenario val- 
idated by simulations. In the ^4^;I case, the last row of 
"heavy" interfacial beads absorbs on a "fast" time scale 
the main part of the energy and linear momentum of the 
incident solitary wave-front (assumed parallel to the in- 
terface) and subsequently decelerates on a "slow" time 
scale, generating a train of (asymptotically) separated 
solitary waves in the "lighter" sonic vacuum. Crucial 
to understanding this phenomenon is the role of contact 
breaking at the interface and the resulting break-down 
of the continuum approximation. When a strongly non- 
linear wave is incident at an oblique angle to the inter- 
face, we find that the angles of refraction and reflection 
are surprisingly well captured by a granular analogue of 
Snell's law that holds irrespective of the solitary wave- 
front amplitude. 

Simulations of parallel interfaces. In order to investi- 
gate solitary waves scattering at two dimensional gran- 
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FIG. 1. Time sequence leading to the generation of a soli- 
tary wave train in simulations. The (red) beads on the left of 
the interface constitute the heavier medium with mass mi 
and the (yellow) beads on the right of the interface con- 
stitute the lighter medium with mass rti2. The mass ratio 
A = 2& = 0.125. The velocity field overlayed in green, de- 
notes the instantaneous speeds of the beads. 



ular interfaces, we performed molecular dynamics simu- 
lations of an impact experiment performed on hexagonal 
packings of 10 4 frictionless spherical grains. As shown in 
Fig. [TJ an interface is introduced by assigning a mass m\ 
to the rows of grains on its left (shown in red) and a mass 
rri2 to rows on its right (shown in yellow). Both portions 
of the hexagonal lattice are comprised of grains with zero 
initial overlap and equal diameters. Two grains of radius 
R and masses {rrii,mj} at positions {x*i,Xj} interact via 
a one-sided non-linear repulsive potential following Hertz 
law GJ 
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(1) 



only for positive comprcssional strains 5y = 2R — \xj — 
0, when 5^ < 0. Here, the inter- 
is expressed in terms of 



xj\ > 0, otherwise Vij 
action parameter 



the effective Young's modulus of the two particles, E*j, 
see Ref. [l3[ for more details. At t = we impart to the 
left-most row a speed u p and subsequently integrate New- 
ton's equations of motion numerically subject to periodic 
boundary conditions perpendicularly to the direction of 
propagation. 

As shown in Fig. [IJa) , this initial condition leads to the 
generation of a non-linear wave front parallel to the inter- 
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face traveling towards the right with a speed V s ~ u p 
analogously to solitary waves in granular chains 0. At 
later times shown in Fig. [ljb), when the solitary wave 
has interacted with the interface, we see a ruptured in- 
terface with one of the interfacial rows of heavy (red) 
beads "dancing" in contact with the lattice of lighter 
(yellow) beads, throttling the generation of an oscillatory 
wave profile in the lighter lattice close to the interface. 
This oscillatory wave is subsequently disintegrated into 
a sequence of separate solitary waves, as shown in Fig. 
[IJc). The separate solitary waves propagate with differ- 
ent speeds (dependent on their amplitude), while a sec- 
ond collision of the "dancing" interfacial row of particles, 
shown in Fig. (Ud), generates a second delayed solitary- 
wave train with smaller amplitudes, sec Fig. QJe) an d 
movie 1 in SI. 

Quasiparticle Collision Model. We take advantage of 
the isotropic elasticity of the hexagonal lattice to assume 
that the dynamics of a solitary wave-front parallel to the 
interface, as in Fig. [TJ is effectively one dimensional and 
governed, in the continuum limit, by the non-linear wave 
equation Q 
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(2) 



where c is a material constant and t) is the strain field 
£{x,t) = —d x u(x,t) expressed in terms of the particle 
displacement field, u(x,t), along the x direction. The 
left-hand side of Eq. @ is the standard inertia term, the 
second term on the right-hand side captures non-linear 
dispersive effects, while the first arises from the restoring 
force as in the wave-equation, if one considers that the 
force is not linear, but it depends on S, 3 ' 2 according to 
Hertz law. 

A strongly nonlinear solitary wave solution of Eq. ([2]) 
can be derived analytically [2j and it has been vali- 
dated by extensive simulations and experiments mostly 
on granular chains @, 0, S Crucially, the total 



energy E 



2m,. 



carried by the solitary wave depends 



l RE ij 



quadratically on its total momentum P, which allows to 
define an effective mass m G ff ~ 1.4m for the solitary wave 

The notion of the solitary wave as a quasiparticle al- 
lows to construct a simple quasi one-dimensional model 
for the generation of the solitary wave train, illustrated 
schematically in Fig. [2] At to, we assume that a chain 
of light yellow beads is uncompressed and all the energy 
and linear momentum P , carried by the incident solitary 
wave, is concentrated in the heavy red interfacial parti- 
cle. At a subsequent time t\ a single solitary wave is 
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FIG. 2. (a) Schematic illustration of our model for the forma- 
tion of a solitary-wave train, side-view, (b) Momentum ratios 
-p^- between the n-th solitary wave and the leading one in 
the train for A — — — 0.125. Red circles are the theoretical 

mi 

predictions while the black squares are the numerical values 
from the simulations of Fig. 1. 



generated in the light chain by reducing the energy and 
linear momentum of the interfacial heavy particle. We 
apply conservation of energy and linear momentum to 
the collision process between the "dancing" bead with 
mass mi and the solitary wave, treated as a quasiparti- 
cle with mass m2, e ff- We calculate the momentum of the 
"dancing" interfacial particle Pi i after the first collision 
as 



Pi,i = 



Po(B-l) 



(3) 



where B = 



(B + l) ' 
The momentum P2,i carried by the 



first leading solitary wave at t = t\ is 

2P 



B + l 



(4) 



At time another independent single solitary wave is 
generated in the "light" chain, further reducing the en- 
ergy and linear momentum of the "dancing" interfacial 
particle. Upon applying conservation of energy and linear 
momentum, as before, and assuming that the first soli- 
tary wave does not participate in this process, we find 
the momentum of the "dancing" particle at t = £2, Pi, 2, 
and of the second solitary wave, P2.2, as 



Pl,2 = 



Po (b - iy 

(B + l) 2 



p2,2 — 



2P (B - 1) 
(B + l) 2 



(5) 



Upon iterating this process n times , we find that the 
"heavy" interfacial bead at t = t„ is left with a linear 
momentum Pi,„ while the n-th solitary wave carries mo- 
mentum P-i^n given by 

Po(P-l)" 2P (P-1)- 1 




time 




FIG. 3. Snapshot of the simulations showing a solitary wave 
incident upon an interface separating two hexagonal lattices 
in a sonic vacuum, (a) For A = — — 0.125, the transmitted 
disturbance propagates in the form of a non-linear oscillatory 
wave analogous to the train of solitary waves shown for gran- 
ular chains in Fig. [T]for A < 1. (b) For A > 1, we find both 
a reflected and transmitted solitary wave as shown in figure 
for A = 3. 



FigpJJb) illustrates the favorable comparison of -pr 1 ^ = 

)n — 1 
against numerical data (red circles) for A = 

0.125. The amplitudes of the delayed secondary sequence 
of solitary waves generated is neglected in our approxi- 
mate model. 

Oblique interfaces The simulation snapshots in Fig. 
[3] illustrate the propagation and interaction of a soli- 
tary wave-front with an oblique interface separating two 
hexagonal granular lattices for the case A = — < 1 in 
panel (a) and for the case A > 1 in panel (b) - see also 
movies 2a and 2b in SI. From these simulations, we have 
determined numerically the angle of refraction (9 re f r for 
different values of the angle of incidence as shown in 
the inset of Fig. 0Ja). 

Inspection of the main panel of Fig. SJa), suggests that 
a linear relationship exists between sin # rG f r and sin f5; for 
mass ratios 4 < 1 (squares) and A < 1 (circles). We 
now show that the measured proportionality coefficient 
is consistent with a non-linear analogue of Sncll's law 
Vosin# re f r = V^.isinf^i where {V2.i,Vb} denote respec- 
tively the speeds of the (leading) refracted solitary wave 
and of the incident one. To work out explicitly the depen- 
dence of on the mass ratio A < 1, we employ a rea- 
soning similar to the one leading to Eq. (QJ that accounts 
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FIG. 4. (a) Angle of refraction # ro f r vs angle of incidence 0i 
for the hexagonal lattice when A < 1. The square (circle) 
symbols correspond to the case when A = 0.125(j4 = 0.9) 
and compares the numerically obtained ratio of the speed of 
the leading transmitted solitary wave to the incident solitary 
wave, against the analytical estimates given by solid curves, 
(b) Comparison of numerical data (symbols) with the analyt- 
ical estimate (solid curves) for the angle of refraction ( black 
data ) and angle of reflection (blue data) vs the angle of in- 
cidence for the hexagonal lattice when A > 1. The insets to 
(a) and (b) , describes the relevant angles, where the interface 
is shown as the dashed (red) line and arrows represent the 
direction of propagation of the solitary wave front (thick dark 
region) . 



for the discrete mechanism at play near the interface. 
Upon making use of the scaling relation P ~ m 3 V s 5 be- 
tween the total momentum P carried by a solitary wave 
and its speed V s j|, we find in the limit of small 0i that 



sin 



sin 



{B + 1)A 3 



(7) 



The right hand-side of Eq. ([7]) is the slope of the con- 
tinuous line plotted for A = 0.125 and A = 0.9 in Fig. 
B|a). It matches the data (open symbols) obtained by nu- 
merically evaluating the left hand-side of Eq. (J7J . This 
agreement corroborates the non-linear analogue of Sncll's 
law. Note, for the case A < 1, the reflected wave ampli- 
tude is negligibly small, see Fig. [3^a). 

By contrast in the case A > 1, Figure [3] (b) shows 
a train of reflected waves with most of the energy con- 



centrated in the leading solitary wave, which allows us 
to neglect rupturing at the interface. As shown in Fig. 
|4] (b) , the angle of reflection (9 ro fl is not equal to the an- 
gle of incidence 8i , since the reflected solitary wave speed 
V^-efl is not the same as the incident speed Vq. Instead, 
we evaluate the ratio in analogy with the derivation 
of Eq. (|7|). Since in this case, the appropriate collision 
model is that between two solitary wave quasiparticles, 
the effective mass contribution cancels out and we obtain 
an equation analogous to Eq. ([3]) with the replacement 
B 1/A and a sign reversal. Thus, Snell's law for re- 
flection assumes the form 



afl) 



sin 



i) 



A-l 



(8) 



The right hand-side of Eq. (jHJ is the slope of the continu- 
ous (blue) line plotted for A = 3 in Fig. HJb) and matches 
the numerical data (blue circles) obtained from evaluat- 
ing the left hand-side of Eq. ([8]). A similar agreement is 
found also for the angle of refraction in this case, as illus- 
trated by the data (black squares) and (black) continuous 
line in Fig. H£b). Despite the fact that our granular lat- 
tices are in a state of sonic vacuum, the ratio between 
the sine of the angles in Snell's law is nearly independent 
of the amplitude of the incident solitary wave front (see 
Figure SI in SI). This seemingly puzzling observation can 
be rationalized by viewing the incident and (leading) re- 
flected or refracted solitary waves as quasiparticles scat- 
tering at the interface whose speeds are proportional to 
the incoming ones. 

Snell's law also implies the existence of a critical an- 
gle of incidence, 9 C , for which the transmitted solitary 
wave will propagate in a direction parallel to the line 
of the interface or # ro f r = 90° in Eq. (JT)), with the re- 
placement B — >• 1/A, In this case, the fronts of the in- 
coming and transmitted solitary wave are no longer con- 
tinuous over the interface. Moreover, the incoming and 
reflected solitary waves cross each other at the interface, 
unlike the case when f9 r efr < 90°, when the two just touch 
each other. Figure S2 in SI shows the delayed reflection 
phenomenon that occurs when the angle of incidence is 
greater than 6 C . 

Conclusion. We have constructed a discrete model 
that predicts the amplitudes of transmitted and reflected 
solitary wave-fronts from a 2D granular interface. We 
find that the angle of refraction and reflection follow an 
analogue of Snell's law in which the solitary wave speed 
replaces the vanishing speed of sound. 
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